Rabi Wave Packets and Peculiarities of Raman Scattering in Carbon Nanotubes, 
Produced by High Energy Ion Beam Modification of Diamond Single Crystals 
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QED-model for multichain coupled qubit system, proposed in was confirmed by Raman scat- 
tering studies of quasi- ID carbon zigzag-shaped nanotubes (CZSNTs), produced by high energy ion 
beam modification of natural diamond single crystals. Multichain coupled qubit system represents 
itself Su-Schriffer-Heeger er-polaron lattice, formed in CZSNTs plus quantized external electromag- 
netic (EM) field. New quantum optics phenomenon - Rabi waves, predicted in Q] has experimentally 
been identified for the first time. It is shown, that Raman spectra in quasi- ID CZSNTs are quite 
different in comparison with well known Raman spectra in 2D those ones. They characterized by 
semiclassical consideration by the only one vibronic mode of Su-Schriffer-Heeger cr-polaron lattice 
instead of longitudinal and transverse optical phonon G + and G~ modes and the out-of- plane radial 
breathing mode, which are observed in Raman spectra of 2D single wall nanotubes. It is consequence 
of 2D - ID transition in all physical properties of nanotubes. It is shown, that strong electron-photon 
coupling takes place in CZSNTs by interaction with EM-field and quantum nature of EM-field has to 
be taken into account. It has been done for the first time in stationary spectroscopy at all. All opti- 
cal spectra, in particular, Raman spectra are registered by usual stationary measurement technique 
in nonequilibrium conditions, which are the consequence of Rabi wave packets' formation. It leads 
in its turn to appearance of additional lines, corresponding to revival part of inversion dependence 
of joint EM-field + matter system in frequency representation. 

PACS numbers: 42.50.Ct, 61.46.Fg, 73.22.-f, 78.67.Ch, 77.90.+k, 76.50.+g 



I. INTRODUCTION 



Quantum electrodynamics (QED) takes on more and 
more significance for its practical application and it, in 
fact, becomes to be working instrument in spectroscopy 
studies and industrial spectroscopy control. QED-model 
for multichain coupled qubit system was proposed in [l[ . 
It is generalization of the model described in , gener- 
alizing, in its turn, Tavis-Cummings model (3j by taking 
into account the ID-coupling between qubits. The most 
substantial result in [2] is the prediction of new quantum 
optics phenomenon - Rabi waves and Rabi Wave packets' 
formation. It is substantial, that in the model, proposed 
in [l[ the interaction of quantized EM-field with multi- 
chain qubit system is considered by taking into account 
both the intrachain and interchain qubit coupling. Tne 
model is presented with an example of perfect carbon 
zigzag-shaped carbon nanotubes (CZSNTs). CZSNTs 
can be considered to be the set of n carbon backbones of 
trans-polyacetylene (t-PA) chains, which are connected 
between themselves. Given rt-chain set can be consid- 
ered to be a single whole, which holds the quasi-one- 
dimensionality of a single chain. It seems to be cor- 
rect for perfect CZSNTs, if their diameter is < 1 nm 
(see Sec. III). At the same time it is well known, that 
free standing nanotubes were considered theoretically to 
be 2£>-strutures and two-dimensional lattice structure 
of a single wall carbon nanotube (SWNT) is specified 
uniquely by the chirality defined by two integers (n, m) 
H3|i 0- Two in-plane G point longitudinal and trans- 
verse optical phonon (LO and TO) modes and the 



out-of-plane radial breathing mode (RBM) [2(| are ob- 
served in the Raman spectra of SWNTs. The LO and TO 
phonon modes at the G point in the two-dimensional Bril- 
louin zone are degenerate in graphite and graphene, how- 
ever they split in SWNTs into two peaks, denoted by G + 
and G~peaks, respectively, [Il[ because of the curvature 
effect. The agreement of experimental Raman studies of 
carbon nanotubes (NTs) with diameter > lnm with 2D 
SWNT-theory unambiguously indicates, that given NTs, 
produced by CVD-methods and like them are really 2D 
systems. At the same time the narrow NTs with diame- 
ter < lnm cannot be considered strongly speaking to be 
2D-systems, they are quasi-lZ? systems, at that 2D-1D 
transition takes place, theoretical explanation of which 
has been done for the first time in 1]. Earlier experi- 
mental results for 2D- ID transition, obtained by electron 
spin resonance (ESR) method in quasi-lD CZSNTs, pro- 
duced by high energy ion beam modification (HEIBM) 
of natural diamond single crystals will be summarized in 
presented work and new experimental results obtained 
by Raman scattering method will be presented. 

The aim of given work is also to represent the experi- 
mental evidence for multichain QED-model, proposed in 
[l| and to confirm the theoretical prediction of Rabi wave 
phenomenon in The new phenomenon predicted is 
quantum coherent effect. It means, that along with re- 
quirement of quasionedimensionality the requirement of 
structural perfectness arises. To observe quantum optical 
coherent effects on NTs, the ensemble of NTs has to be 
homogeneous. It means, that any dispersion in axis di- 
rection, chirality, length and especially in diameter both 
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for single NT along its axis and between different NTs 
in ensemble has to be absent, axial symmetry has to be 
also retained, that is, there are additional requirements 
in comparison with, for example, t-PA technology. The 
CVD-technology of NTs production and many similar to 
its ones seem to be not satisfying to above-listed require- 
ments at present. It means, that experimental results and 
their theoretical treatment will be different in both the 
cases, that really takes place (see for more details Sec. Ill 
and Sec. IV). The situation seems to be analogous to some 
extent to the solid state physics of the same substance 
in single crystal and amorphous forms. The technology, 
based on HEIBM, satisfy given requirements. 

The paper is organized in the following way. In Section 
II, experimental technique and Raman scattering results 
are described. In Section III, the comment and some 
development of Su, Schrieffer, Heeger (SSH) model for 
organic conductors are given. In Section IV, experimen- 
tal results are discussed. In Section V, the conclusions 
are presented. 



II. EXPERIMENTAL TECHNIQUE AND 
RESULTS 

Samples of type Ha natural diamond, implanted by 
high energy ions of copper (63 MeV, 5 x 10 14 cm -2 ) 
and boron (13, 6 MeV) have been studied. Ion implanta- 
tion was performed along (111) crystal direction. Raman 
scattering (RS) spectra were registered in backscattering 
geometry. Laser excitation wave length was 488 nm, rect- 
angular slit 350x350(/im) 2 was used, scan velocity was 
100 cm" 1 pro minute. 

The spectra observed are presented in Figures 1 to 3. 
All the spectra are 180 out of phase, that is dip positions 
correspond to maxima of the signals. The spectrum, 
presented in Figure 1, is characteristic of the only ion 
beam modified region (IBMR) of the sample, since the 
RS-line near 1332 cm' 1 , characteristic for diamond sin- 
gle crystals, is absent in the spectrum. The RS-lines with 
peak positions 656.8±0.2 cm" 1 , 1215±1 cm" 1 , 1779.5±1 
cm' 1 and 2022. 3±0. 5 cm' 1 correspond to optical cen- 
ters in IBMR, which are active by laser excitation trans- 
versely to sample surface from implanted side. At the 
same time, by laser excitation of the same sample from 
opposite unimplanted side the qualitatively other picture 
is observed, see Figure 2. Firstly, the line, characteristic 
for diamond single crystals, is now presented in the spec- 
trum (the whole amplitude of given line is not shown 
in Figure 2). It is interesting, that its frequency value 
(1328.7 cm' 1 ) is slightly shifted from usually observed 
value near 1332 cm' 1 . It indicates on renormalization of 
optical phonon by Cu-implantation in the whole diamond 
matrix. 

Secondly, RS-lines with peak positions at 354.6, 641.8, 
977.1 (±1 cm' 1 ), 1274.1 ±2 cm' 1 and more weak pro- 
nounced lines at 1569 ±3 cm' 1 , 1757±5 cm' 1 were ob- 
served. Moreover, all the lines listed are superimposed 
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FIG. 1: Spectral distribution of Raman scattering intensity 
in diamond single crystal, implanted by high energy copper 
ions, the excitation is from implanted side of the sample. 

now with very broad (its linewidth is 1720±20 cm' 1 ) 
asymmetric line with peak position 116Q±10em . 

The spectrum of boron implanted sample, Figure 3, 
was measured in the range (1000 - 2100) cm' 1 . It is seen, 
that the spectra of boron and copper implanted samples 
are qualitatively similar in given spectral range. However 
numerical values - 1212.3±1, 1772. 5±1, 2011±0.5cm" 1 
- of the peak positions are slightly different, they are 
shifted to low frequency region. It is interesting, that 
there is regularity in peak position shift, which is increas- 
ing with frequency increase, and it is equal to 2.8, 7, 11.3 
cm' 1 correspondingly. 

It is strong indication, that all three lines belong to the 
same optical system in both the samples studied. The 
line with peak position 1331.95 zbO.lcm" 1 was also pre- 
senting. Its frequency value coincides with the value of 
optical phonon peak in conventional natural diamonds. 
The presence of given RS-line by excitation from im- 
planted side can be determined by two factors. Firstly, 
the effective thickness of IBMR is substantially less, it 
is ~ 1.4/zto, in comparison with Cu-implanted sample 
(~ 7/im) [Ty] ! consequently, laser excitation can reach an 
unimplanted region. It can also be suggested, that the 
modification is not entire in near surface region. Given 
results seem to be the first results in application of RS- 
spectroscopy to study of IBMR in diamonds. Let us also 
to pay attention, that the length of CZSNTs is deter- 
mined by the implantation with the ion energy ~ lMev 
per nucleon by the thickness of IBMR, that is it gets 
in usual for free standing NTs micron-th range. It fol- 
lows from ESR experiments on layer- by- layer removal of 
IBMR, see [ll[ and the references therein. 

Comparing RS-spectra, represented in Figures 1 to 
3 with well known spectra of NTs produced by CVD- 
methods it is seen, that difference between the optical 
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FIG. 2: Spectral distribution of Raman scattering intensity 
in diamond single crystal, implanted by high energy copper 
ions, the excitation is from unimplanted side of the sample. 
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FIG. 3: Spectral distribution of Raman scattering intensity in 
diamond single crystal, implanted by high energy boron ions, 
the excitation is from implanted side of the sample. 



characteristics of the NTs, produced by HEIBM of dia- 
mond single crystals and between the NTs, produced by 
other methods, including methods of ultra-small NTs' 
production in zeolite matrix, is greatly, while there is the 
essential similarity of ESR-data with those ones for ultra- 
small NTs, embedded in zeolite matrix (see Sections III, 
IV). 

On the other hand, comparing the figures 1 and 2 be- 
tween themselves, we can straight away conclude, that 
RS-lines observed cannot be attributed to usual atomic or 



molecular vibrations including the local vibration modes 
like to well known G-bands and RBM, observed in NTs, 
produced by CVD-methods, since the RS-spectra, corre- 
sponding to usual local vibrations cannot be dependent 
on the change of propagation direction of exciting light 
into opposite. 

Therefore 2D-theory, developed for explanation of op- 
tical properties of 2D-NTs cannot be applicable for per- 
fect NTs of small diameter. It is strong argument for 
application of the QED-model, proposed in [l| for cor- 
rect description of the objects studied. 

Since the theoretical model of small quasi-lD CZSNTs 
[U is based partly on SSH-model of ID organic conju- 
gated conductors let us consider SSH-model more care- 
fully. It will be shown, that it needs in some corrections 
and in development. 

III. COMMENT TO SSH-MODEL FOR 
ORGANIC CONDUCTORS AND SOME ITS 
DEVELOPMENT 

It is substantial, that the physical properties of t-PA 
chains are studied very well both theoretically and exper- 
imentally. Let us touch ongiven subject in more detail. 
Su, Schrieffer, Heeger [J have found the most simple 
way to describe mathematically the chain of t-PA by con- 
sidering it to be Fermi liquid. It is now well known SSH- 
model. The most substantial suggestion in SSH-model is 
concerned the number of degrees of freedom. Su, Schri- 
effer, Heeger suggested, that the only dimerization coor- 
dinate u n of the n-th Cff-group, n = 1, N along chain 
molecular-symmetry axis x is essential for determination 
of main physical properties of t-PA. Other five degrees of 
freedom were not taken into consideration. Nevertheless, 
the model has obtained magnificent experimental confir- 
mation, that requires additional argumentation. At the 
same time it will be shown in the next subsection, that 
SSH-model itself needs in some corrections, which being 
to be of the principle lead nevertheless qualitatively to 
analogous conclusions, concerning physical properties of 
the system studied. 

A. Corrections to SSH-Model 

We have repeated the calculation, presented in [3| , Q . 
It has been found, that there is the second solution, which 
leads to independent branch of quasiparticles. Let us 
show it. We preserve all designations for quantities from 
5]- 

The expression for Hamiltonian in reduced zone has 
been obtained in the following form 

= EE^* " + ^ k k A k }(n c ks - nls) 
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where ak and (3k are determined by 



a. The first condition 
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are occupation number operators in C- and 
1^-bands correspondingly, e is unit operator. The known 
SSH-solution corresponds to choice of lower signs in rela- 
tions ©. The choice of upper signs gives independent so- 
lution. Therefore, for the energy of quasiparticles E^ u \ 

E k ^ in C- and F-bands, corresponding to upper-sign- 
solution we obtain 



E, 



^1 



Sn 



v[u] _ 6EM(u) 
k " ML 



Ek 

~ ( E k h 



(5) 



We see, that it differs from the energy of quasiparticles 
E^ l \ E% in C- and TA-bands, corresponding to lower- 
sign-solution (SSH-solution), which is 
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£M(u) and £?M( U ) in © and © are eigenvalues of op- 
erator H(u), which correspond to upper and lower signs 
in @, n ks , n v ks are eigenvalues of operators of particle 
numbers in C-band and V^-band correspondingly. The 
coefficients ak and (3k were found in Q from the condi- 
tions of energy minimum. However, the only necessary 
conditions were used. At the same time the sufficient 
conditions for the minimum are substantial in given case, 
they change the role of both solutions. Really, since ak 
and (3k are coupled by the condition 
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for the SSH-solution and 



£fc 
E k 


s A l 
< ~Ek~ 


\{n c ks - 




£fc 
E k 


A 2 


K« - 


-n\ s ) 



(8) 



{£ fe (l + ^<^|K fes -nD<0}, 

&k &k 



(9) 



for the additional solution. It is seen, that the first condi- 
tion is realizable for the quasiparticles of both the kinds, 



< 0) and in nonequilib- 



at that in equilibrium (n ks — n v ks 
rium (n ks — n v ks > conditions. 

b. The second condition The second condition is the 
same for both the solutions and it is 
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c. The third condition For the SSH-solution we have 
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It means, that SSH-solution is unapplicable for de- 
scription of standard processes, passing near equilibrium 
state by any parameters. The quasiparticles, described 
by SSH-solution, can be created the only in strongly 
nonequilibrium state with inverse population of the levels 
in C- and y-bands. At the same time for the solution, 
which corresponds to upper signs in ([2]), we obtain 
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which is realizable both in near equilibrium and in 
strongly nonequilibrium states of the systems by corre- 
sponding choose of parameters. 

Let us consider the continuum limit for the ground 
state of the i-PA chain with quasiparticles of given 
branch. Taking into account, that in ground state n ks = 
0, n ks = 1 we have 



we have conditional extremum and sufficient conditions 
for the minimum can be obtained by standard study. The 
second differential of the energy to be the function of 
three variables ak, (3k and A, where A is coefficient in 
Lagrange function for conditional extremum, has to be 
positively defined quadratic form. From the condition of 
positiveness of three principal minors of quadratic form 
coefficients we obtain the following three sufficient con- 
ditions for the energy minimum 
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then, calculating the integral, we obtain 
4Nt 



1 + z 2 



'-{F(-,l-z 2 ) + 

IT 2 



— -[£(|, 1 - z 2 ) - F(|, 1 - z 2 )]} + 2NKu\ 



(13) 



(14) 



5 



where -F(§, 1 — z 2 ) is the complete elliptic integral of the 
first kind, E(^,l — z 2 ) is the complete elliptic integral 
of the second kind, z 2 = -^p. Approximation of (H"4l at 
z<l gives 
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It is seen from (|15p . that the energy of quasiparticles, 
described by solution, which corresponds to upper signs 
in ([2} has the form of Coleman- Weinberg potential with 
two minima at the values of dimerization coordinate uo 
and —uo like to energy of quasiparticles, described by 
SSH-solution @. 

Therefore, all qualitative conclusions of the model pro- 
posed in [j| are holding, however for the quasiparticles, 
corresponding to the second-branch-solution. 



B. Slater Principle and SSH-Model 

In [l[ was suggested, that success of SSH-model is the 
consequence of some general principle and it was shown, 
that given general principle is really exists. Main idea 
was proposed by Slater at the earliest stage of quantum 
physics era already in 1924, that is before creation of 
quantum mechanics and quantum electrodynamics. It 
is - "Any atom may in fact be supposed to communi- 
cate with other atoms all the time it is in stationary 
state, by means of virtual field of radiation, originating 
from oscillators having the frequencies of possible quan- 
tum transitions..." @. The development of given idea is 
based on the results of the work [7]. It has been found 
in 0], that Coulomb field in ID-systems or 2D-systems 
can be quantized, that is, it has the character of radia- 
tion field and it can exist without the sources, which have 
created given field. Consequently, Coulomb field can be 
considered to be "virtual" field in Slater principle and it 
can be applied to both t-PA and to quasi-lD-NTs. It 
produces in t-PA the preferential direction in atom com- 
munication the only along chain axis (to be consequence 
of quasi-one-dimensionality). It is reasonable to suggest, 
that given direction remains to be also preferential by in- 
teraction with external EM-field, then the explanation of 
the success of SSH-model, taking into consideration the 
only one degree of freedom of n-th Cff -group instead of 
six ones becomes to be natural. It determines also the 
applicability of SSH-model in the part concerning active 
degrees of freedom by interaction with external EM-field 
to quasi-one-dimensional CZSNTs, while the interaction 
between the chains, which produce CZSNTs is strong (see 
further). 

Moreover Slater principle can be applied to 2D- 
SWNTs. However the physical consequences of Slater 
principle application are quite different for quasi- ID and 
2D systems. In quasi-lD-systems, that is in t-PA and in 



quasi-lD SWNTs Coulomb field can be considered to be 
"virtual" field with propagation direction the only along 
t-PA chain and NT-axis correspondingly. In other words 
it produces preferential direction in atom communication 
the only in one direction (to be consequence of quasi-one- 
dimensionality) , and given direction remains to be pref- 
erential by interaction with external EM-field, that ex- 
plains qualitatively the success of SSH-model in the sense 
that degrees of freedom, realized by bonds, projections of 
which are not coinciding with chain molecular-symmetry 
axis direction, can really be not taken into considera- 
tion for experiments with the participation of external 
EM-field. The consequence of given conclusion is the 
following. Quasi-lD CZSNT can be represented by the 
model consisting of the set of n interacting between them- 
seves equivalent carbon backbones of t-PA chains, while 
in real quasi- ID CZSNTs the adjacent chains are rep- 
resenting the mirror image of each other (relatively the 
corresponding radial plane). Given detail of SSH-model 
is used further in the model of quasi-lD CZSNTs. At the 
same time, there are existing in 2D-systems including 2D- 
SWNTs two preferential directions of Coulomb "virtual" 
field propagation. It means, that for 2D-SWNTs all de- 
grees of freedom, which are relevant to the bonds in rolled 
graphene sheet will be essential by interaction with ex- 
ternal EM-field and the simplification, which is fruitfully 
used in 1D-SSH model becomes to be incorrect for 2D- 
SWNTs in full accordance with existing 2D-SWNTs the- 
ory and Raman scattering experiments, see for instance 

m Q, M, M- 

Experimental confirmation for applicability of Slater 
principle to quasi-lD-CZSNTs and for the model of 
quasi-lD-CZSNTs above proposed follows from ESR- 
studies in rather perfect quasi-lD-CZSNTs, produced by 
HEIBM of diamond single crysatals, by which the ap- 
pearance of Peierls transition and neutral paramagnetic 
(with spin S = 1/2) SSH-7r-soliton formation were es- 
tablished [Tl|, [l2[. By the way, it means, that perfect 
quasi-lD CZSNTs, characterized by (m, 0) indices, will 
have bandgap like to classical semiconductors at any m, 
including the case m = 3q, q E N, for which the 2D- 
theory, existing at present, predicts the metallic proper- 
ties. It follows also from Slater principle the following. 
Longitudinal and transverse optical phonon graphite-like 
G-modes [jjj], undergoing in 2D-SWNTs splitting into 
G + and G~ modes, respectively, [2l| because of the cur- 
vature effect and the out-of-plane RBM [20(, that is all 
the modes, which are observed in Raman spectra of 2D- 
SWNTs, have to disappear in perfect quasi-lD CZSNTs. 
Really, the fact, that neutral SSH-7r-solitons are respon- 
sible for ESR-spectra in quasi-lD CZSNTs incorporated 
in diamond matrix is direct confirmation to given conclu- 
sion. Neutral (zero charged) SSH-7r-solitons are optically 
inactive [29j . Consequently all 7r-subsystem will be inac- 
tive in optical spectra of quasi-lD CZSNTs. It seems to 
be one of the most substantial characteristics of 2D-1D 
transition in physical properties of CZSNTs. 

It seems to be interesting to predict, what kind of lines 
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have to appear in Raman spectra instead of G + , G~ and 
i?£?-modes by 2D-1D transition. The model of quasi-lD 
CZSNTs proposed in [l[, experimental results on ESR- 
studics of quasi-lD CZSNTs and on optical studies of 
related carbon chain material - carbynes allow to obtain 
given prediction a priori without detailed analytical cal- 
culation. It is sufficient to take into consideration, that 
SSH-model along with the physical basis of the existence 
of solitons, polarons, breathers, formed in 7r-electronic 
subsystem (Tr-solitons, 7r-polarons, 7r-breathers) contains 
in implicit form also the basis for the existence of simi- 
lar quasiparticles in er-electronic subsystem, that is SSH- 
model can be developed. It was done in [27| and in [351 ] - 



C. CT-Quasiparticles in SSH-Model 

The origin of quasiparticles' formation in er-electronic 
subsystem is the same two-fold degeneration of ground 
state of the whole electronic system, energy of which 
in ground state has in correspondence with results in 
the first subsection the form of Coleman- Weinberg po- 
tential with two minima at the values of dimerization 
coordinate u and — u for both kinds of quasiparticles 
in SSH-model. Really the appearance of uq ^ and 
— uq ^ indicates on the alternation in interatomic dis- 
tance. It means, that simultaneously with 7r-subsystem, 
cr-subsystem will also be dimerized. 

The shapes, for instance, of 7r-solitons and tr-solitons 
can be given by the expression with the same mathemat- 
ical form 
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^ (a 



_ sech2! (n-n )a 
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«7r(o-)*J COS — , (16) 



where n, rig are variable and fixed numbers of CiJ-unit in 
C-ff-chain, a is C — C interatomic spacing projection on 
chain direction, v w i a \ is 7r(cr)-soliton velocity, t is time, 
£„.( CT ) is 7r(er) coherence length. It is seen, that 7r-solitons 
and cr-solitons differ in fact the only by numerical value 
of coherence length. Given difference can be evaluated 
even without numerical calculation of the relation, which 
determines the shift of ground state energy of extended 
system by presence of localized perturbation. Actually it 
is sufficient to take into account the known value of £ w 
and relationships [33| 
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where Aoer, Ao,r are a— and 7r-bandgap values at T = 
OK, vp is Fermi velocity. Theoretical value in t-PA 
is 7a, and it is low boundary in the range 7a — 11a, 
obtained for £„• from experiments [2Qj j . Taking into ac- 
count the relationships (|17[) . using the value s» 8.8, 
which was evaluated from t-PA band structure calcula- 
tion in 31], and mean experimental value of coherence 
length £„. = 9a we obtain the value £ CT w 0.125nm. It 
means, that halfwidth of space region, occupied by a- 
soliton in t-PA is « 0, 25nm, that is SSH-cr-solitons are 



much more localized in comparison with SSH-7r-solitons. 
Similar conclusion takes place for SSH-tr-polarons repre- 
senting itself the soliton-antisoliton pair. SSH-c-polarons 
have recently been experimentally detected in the work 
27J , where the formation of polaron lattice (PL) was es- 
tablished. It was found, that two components of each 
elementary unit, that is, of each polaron, possess by two 
equal in values electical own dipole moments, propor- 
tional to spin, which was called electical spin moments 
(ESM), with opposite directions. It was shown, that ex- 
perimental results agree well with PL-formation, which 
means in fact the formation of antiferroelectrically or- 
dered lattice of quasiparticles. Given lattice consists of 2 
sublattices, corresponding to soliton and antisoliton com- 
ponents of polaron. Corresponding chain state is opti- 
cally active and it is characterized by the set of lines in 
infrared (IR) spectra, which were assigned with new op- 
tical phenomenon - antiferroelectric spin wave resonance 
(AFESWR). Let us remember that carbynes are organic 
quasi-one-dimensional conductors with the simplest, con- 
sisting the only of the carbon atoms, chain structure. At 
the same time the presence of two electronic tt x and ir y - 
subsystems, which are "hung" on a single cr-subsystem 
means, that the ground electronic state is similar to two- 
dimensional Coleman- Weinberg potential with four min- 
ima at the values of dimerization coordinate uo and —uq. 
In other words, ground electronic state in carbynes is 
four-fold degenerate, which leads to a substantially more 
rich spectrum of possible quasiparticles, discussed in (36j 
and in [28| . 

Taking into accout that cr-subsystems are very similar 
in carbynes and in t-PA and consequently in an arbitrary 
chain of n-chained quasi-lD CZSNTs, we can evaluate 
the numerical characteristics of AFESWR in in t-PA and 
in quasi-lD CZSNTs. Really, central mode, that is AFR 
mode has the frequency value v°{C) in carbyne sample 
equaled to 477 cm" 1 , the splitting parameter in IR de- 
tected AFESWR spectra was equal to 150 cm" 1 . Given 
values of v°(C) and AFESWR-splitting parameter in car- 
bynes allow to estimate the range for expected values of 
analogous parameters in t-PA and in quasi- ID CZSNTs 
in the following way. AFESWR-splitting parameter is 
determined by exchange integrals in cr-electronic subsys- 
tem [27], which seems to be practically the same in car- 
bynes and in t-PA and in quasi-lD CZSNTs, since the 
role of quite different 7r-subsystems in carbynes and in 
t-PA and quasi-lD CZSNTs can be neglected to a first 
approximation. Consequently, the value of AFESWR- 
splitting parameter in quasi-lD CZSNTs and in t-PA has 
to be close to 150 cm -1 and to 300 cm -1 by IR- and RS- 
AFESWR-detection correspondingly 30]. The frequen- 
cies z/£(C) and i£(t - PA), v°{NT) of main AFESWR- 
mode in SSH-a-polaron lattice in carbyne and in t-PA 
and in quasi-lD CZSNTs depend on intracrystalline field 
27], that means, that their values will be different. How- 
ever the values of v°(t—PA), v°{NT) can be evaluated, if 
to take into account the known relation for the vibration 
frequencies of similar centers, the fact of cr-polaron and n- 
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soliton lattice formation in carbynes, leading to change in 
effective masses [28| , presence of two 7r-subsystems in car- 
bynes, difference of coherence lengths in accordance with 
@ of cr-solitons (£„) and 7r-solitons (£,r), band structure 
data for t-PA [3l| and carbynes [32J. 

We have obtained the following frequency ranges for 
IR SSH-cr-polaron lines in t-PA and in quasi-lD CZSNTs 
v a v {t-PA) G (386.7, 603) cmT 1 and v° (NT) g (402.5, 
627.6) cmT 1 . Known IR-mode with the frequency near 
540 cm,- 1 , in t-PA [H| gets to interval (386.7, 603) cm -1 
and can represent itself AFR mode in cr-polaron lattice, 
that is, there is alternative interpretation of given IR- 
mode, ascribing earlier to Goldstone SSH-7r-soliton vi- 
bration mode 29]. It follows from results of [34|, that 
the same spectral interval with slightly different right- 
hand value, equaled to 673.7 cm" 1 , is the evaluation for 
the frequency of cr-polaron main AFESWR-modes in t- 
PA and in quasi-lD CZSNTs, which are RS-active. The 
calculation in |34| does not take into consideration the 
soliton and polaron formation. However cr-polaron for- 
mation does not violate the symmetry of task, that allows 
to conclude, that the same asymmetry value will be re- 
tained in IR and RS spectral distributions. 

Therefore the qualitative semiclassical consideration 
predicts a priori, that Raman spectrum of quasi-lD CZS- 
NTs will consist of the only one line, representing itself 
AFR-modc of cr-polaron lattice with peak frequency po- 
sition in the range (386.7, 673.7) cm -1 . It can be splitted 
into series of AFESWR-modes with average splitting pa- 
rameter near 300 cm -1 . The possibility of AFESWR- 
splitting depend on experimental geometry conditions 
(that is, whether is an antiferroelectric spin wave modes' 
excitation allowed by the experimental geometry or not). 

It will be shown further, that the prediction has to be 
completed, if to take into account the quantum nature of 
EM-held. 

Experimental confirmation for conclusion on 2D-1D 
transition and its main physical properties predicted fol- 
lows from ESR-studies, performed earlier in rather per- 
fect CZSNTs, produced by HEIBM of diamond single 
crysatals [13] , [H| , [13 j and from Raman scattering data, 
presented in Section II, obtained on the same samples. 
Let us reproduce the formulation of the model of quasi- 
lD CZSNTs, given in Q. Quasi-lD CZSNT represents 
itself autonomous dynamical system with discrete circu- 
lar symmetry consisting of finite number n G N of car- 
bon backbones of t-PA chains, which are placed period- 
ically along transverse angle coordinate. Longitudinal 
axes {xi},i = l,n, of individual chains can be directed 
both along element of cylinder and along generatrix of 
any other smooth figure with axial symmetry. It is taken 
into account, that Slater principle like to SSH-model for 
t-PA allows to consider to be active the only degree of 
freedom along axes {xi},i — l,n, of individual chains 
or, in other words, along single hypercomplex axis in a 
hypercomplex number language of the task formulation. 
It is the reason, that the adjacent chains, which repre- 
sent themselves a mirror of each other in real structure, 



becomes to be equivalent in the model, that is, they will 
be indistinguishable, since the only one degree of free- 
dom - dimerization coordinate u m of the m-th C-atom, 
m = l,N along chain molecular-symmetry axis x is sub- 
stantial for determination of main physical properties in 
the frames of the model proposed. 



IV. DISCUSSION 

A. Summarization of Results on HEIBM-Method 
of Nanotubes Production 

Let us give short review concerning HEIBM-mcthod of 
production of incorporated carbon NTs and to represent 
the exact experimental proof of their real formation in 
diamond matrix. 

The first report on the discovery of new carbon phase 
- carbon nanotubes, incorporated in diamond matrix, is 
related to 1990, and it was made during 1990 IBMM- 
Confcrence, Knoxwille, USA. Similar report was also 
represented at E-MRS 1990 Fall Meeting, Strasbourg, 
France, @, that is, the conclusion on new carbon phase 
formation became to be known substantially prior of now 
well known japan discovery of free standing NTs, re- 
ported in [37j, which is related to 1991. Basis experi- 
mental method was ESR. It was found, that central place 
takes in ESR spectra of diamond single crystals, modi- 
fied by high energy ion implantation, the single inten- 
sive line with very unusual radiospectroskopic properties. 
It was anisotropic, however the anisotropy was weak in 
comparison with the anisotropy of point centers in dia- 
mond. Especially interesting, that along with g- value the 
linewidth AH PP was also found to be tensor quantitity, at 
that <?-tensor and A£f pp -tensor were characterized by the 
same axial symmetry group (To| . Moreover the only one 
eqivalent configuration in diamond lattice was presented 
with g\ 1 , AH PP principal directions of axial g- and AH PP - 
tensors, which were strictly coinciding with ion beam di- 
rection, at that <7|| is minimal g- value, AH PP is maximal 
AH pp -value in g and AH PP angular dependences. Other 
lattice equivalent configurations were absent. The kind of 
symmetry group, corresponding to g-tensor and AH PP - 
tensor symmetry was strongly dependent on the choice 
of ion beam direction relatively the crystallographic axes 
of diamond lattice. It was shown for the first time in 
radiospectroscopy, that g-tensor and A£f pp -tensor sym- 
metry corresponding to the most intensive line in the 
spectrum observed are the mapping of the symmetry of 
lengthy objects with macrosizes along implantation direc- 
tion [10| . Detailed studies of ESR spectral angular de- 
pendences in appropriate crystallographic planes allowed 
to establish, that the structures formed by (lll)-HEIBM 
in diamond have tracklike cylindrical shape, which are 
lengthy strictly in (lll)-direction, with the size varying 
from several mkm to several tens of mkm, depending on 
ion energy used. In other words, the first studies in 1990 
allowed to establish the formation of carbon nanotubes 
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of cylindical symmetry shape incorporated in diamond 
matrix. Simultaneously the formation of quite differ- 
ent nanotubes was found. They represent themselves 
crimped cylinder with crimping corresponding to four- 
petal structure in the cross section. They are produced 
by the implantation direction, coinciding with (100) axis 
of diamond lattice 10]. 

The next step was the additional confirmation of the 
structure of cylindical symmetry shape to be really (111)- 
incorporated nanotubes with rolled up graphene sheet in 
zigzag shaped configuration (but not for example cylindi- 
cal rod). It was done by the study of radiospectroscopic 
properties of spin curriers in given tubes. It was found, 
that the system of paramagnetic centers (PC), which are 
responsible for appearance of strong single line absorp- 
tion above described is non-Blochian system. A number 
of distinctive peculiarities have been observed for the first 
time in radiospectroscopy at all. The main ones among 
them are the following. 

1 .It was found, that the shape of resonance lines cannot 
be described by Lorentzian, Gaussian, or by their convo- 
lution. The shape is characterized by essentially more 
slow decrease of the absorption intensity on the wings in 
comparison with known shape functions. It was called 
super-Lorentzian. It was proved, that super-Lorentzian 
is genuine lineshape of resonance absorption (that is, it 
does not represent the superposition of Lorentzians with 
different linewidths) (Tlj . 

2. It was established, that PC-system is nonsaturating 
with superlinear absorption kinetics (that is, the depen- 
dence of resonance signal amplitude on the amplitude of 
magnetic component of microwave field is superlinear). 
Especially interesting, that superlinear absorption kinet- 
ics has been observed by both the registration of res- 
onance signal in phase with high-frequency modulation 
and in quadrature with high-frequency modulation [Io| . 

3. It was found unusual for nonsaturating resonance 
systems dependence of resonance signal amplitude on the 
value of modulation frequency, which was characterized 
by substantial increase (instead decrease) of resonance 
signal amplitude with modulation frequency increasing 

Ei, EH- 

4. It was established the effect of appearance of "phase 
angle" , characterizing the absorption process. It consist 
in that, that although absorption kinetics is nonsaturat- 
ing, the maximum of absorption is achieved not strictly 
in phase with modulation field, that is at modulation 
phase value, which is nonequal to zero. "Phase angle" 
is anisotropic in general case and maximum of its value 
was found to be equal 20 fill ] . 

The analysis of given peculiarities [ll|, [l2[ has led 
to conclusion, that PC, which are responsible for strong 
non-Blochian absorption, are mobile and are character- 
ized by two times of spin-lattice relaxation, by T\ , which 
is comparable with relaxation time of usual point PC 
in diamond single crystals (10 -3 — 10 _4 s) and by very 
short time r (upto 10 _13 s) of conversion of the energy 



of the spin system into mechanical kinetic energy of PC 
motion, which can be considered to be consequence of 
fundamental in magnetic resonance phenomenon of gy- 
romagnetic coupling between magnetic and mechanical 
moments. Therefore T-process is very fast process in 
comparison with any process, realized by means of pure 
magnetic interaction, the relaxation times of which are 
not shorter, than 10 _10 s The possibility of given 
conversion can be realized naturally, if the energy of the 
motion activation is very small, it is typical for mobile 
topological solitons like them ones, identified in t-PA. So 
the conclusion, that the strong non-Blochian absorption 
is determined by mobile topological solitons was obtained 
by natural way. All subsecuent studies have confirmed 
given conclusion. It was found, that there is numerical 
coincidence of the characteristics of SSH topological ir- 
solitons in t-PA and in (111) incorporated NTs. So, the 
values of g-tensor components in Cu-implanted sample 
are g\ — 2.00255 (it is minimal g- value and it is g\\ prin- 
cipal direction), gi = 93 = g± = 2.00273, the accuracy 
of relative g- value measurements is ±0.00002 [l(|. g- 
value of paramagnetic 7r-solitons in trans-polyacetylene, 
equaled to 2.00263 13|, gets in the middle of given rather 
narrow interval of g- value variation of PC in ion produced 
NT. Although anisotropy of paramagnetic 7r-solitons in t- 
PA, which are also considered to be mobile PC, mapping 
the distribution of 7r-electron density along whole indi- 
vidual t-PA chain, is not resolved by ESR measurements 
directly (which in fact is the indication, that chemically 
produced t-PA is less perfect in comparison with NTs in 
diamond matrix), there are indirect evidences on axial 
symmetry of ESR absorption spectra in t-PA too, 
fig . Consequently, the value 2.00263 is mean value and it 
coincides with accuracy 0.00002 with mean value of afore- 
cited principal g-tensor values of PC in NTs. Given coin- 
cidence becomes to be understandable now, if to take into 
account the results of [l|, where is shown, that model of 
CZSNTs, produced by HEIBM of diamond single crystals 
is very similar to SSH-model of t-PA. The main differ- 
ence of both the models consist in the following. Carbon 
chain backbone of t-PA single chain can be considered to 
be ID-object in the space 3D <g> Zi, carbon chain back- 
bone of single CZSNT can be considered to be ID-object 
in 3D ® Z n space. Here 3D is real Euclidian space, Z n is 
hypercomplex commutative ring, which is direct sum of 
n complex spaces C. 

z n = c (sc e ...®c, (is) 

Z\ = C. We see, that average ESR properties of ID- 
objects in 3D ® Z\ space and in 3D ® Z n space are iden- 
tical with high precision. 

On the other hand we see, that it was possible to de- 
termine anisotropy of g-values with very high precision, 
that indicates on the very perfect and homogeneous axi- 
ally symmetric NTs with symmetry axis strictly along ion 
beam direction and to establish the origin of PC in ion 
produced NTs to be paramagnetic topological 7r-solitons 
(SSH-solitons), considered to be mobile PC, mapping the 
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distribution of 7T- electron density along whole individual 
NT [ll|, [Hj]. Thereby immediately from ESR studies 
was found, that structural element of new carbon phase 
with cylindrical symmetry produced in diamond by (111) 
high energy implantation is t-PA chain backbone. At 
the same time t-PA chain backbone is basic element for 
graphene, graphite, NTs formation and also for the car- 
byne formation. Carbynes were studied by ESR, they 
characterized by quite different ESR-spectra [26J. The 
possibility of formation of cylindrical graphite rod is also 
not corresponding, since graphite is characterized by the 
other ESR spectra. So it was concluded unambiguously, 
that the spectra observed by (111) high energy implanta- 
tion are corresponding the only to graphene sheet, which 
is rolled up in that way, in order to t-PA chain backbones 
were directed along (lll)-axis of diamond lattice, coin- 
ciding with ion beam direction, that results in CZSNTs 
production. 

Therefore the task of identification of electronic and ge- 
ometrical structure of CZSNTs was solved by ESR spec- 
troscopy methods, which are the most powerfull methods 
[it seems to be commonly accepted opinion] in the struc- 
tural studies of nanosize objects. Thousands spectra, 
including the spectra reported in literature for natural 
diamonds single crystals both in starting materials and 
irradiated with electrons, neutrons, ions and all known 
ESR data in related carbon materials - polycrystalline 
diamond films, fullerenes, coals and graphite were stud- 
ied and analyzed. 

It is substantial, that the CZSNTs are not random in 
orientation. They are produced strictly in ion beam di- 
rection and they have usual length by ~ lMev per nu- 
cleon ion energy, it is proved in details in the works above 
cited. Naturally random processes are also presenting 
by HEIBM. The knocked diamond atoms produce point 
defects in lateral to ion track directions, their concen- 
tration distribution in the sample represents itself the 
overlap of random distributions along directions, which 
are perpendicular to each single ion track axis. They 
were also detected and studied by ESR, see, for instance 
[T3 |. The main difference of our results on NTs from 
well known ones is based on other dimensionality of CZS- 
NTs studied. We are dealing with quasi-one-dimensional 
CZSNTs which have to possess and really possess by 
quite other physical properties including optical proper- 
ties in comparison with two-dimensional CZSNTs, which 
have naturally graphite-like optical spectra. In addition 
we see also, that the spectra observed by Raman spec- 
troscopy in quasi-one-dimensional CZSNTs can be pre- 
dicted qualitatively, if to take into consideration the only 
ESR data (optical data in carbynes were used addition- 
ally for quantitative evaluation of optical characteristics, 
that is AFESWR characteristics). 

It was above indicated, that free standing nanotubes 
were considered theoretically the only to be 2Z?-strutures. 
At the same time there are along with experimental 
ESR results on SWCNTs incorporated in diamond ma- 
trix above summarized, the experimental ESR and Ra- 



man scattering results on narrow free standing tubes with 
diameter < lnm, indicating , that they cannot be consid- 
ered strongly speaking to be 2D-systems, they are quasi- 
1D systems. Really, in [22j is reported on the develop- 
ment of SWCNTs of 0.4 nm diameter - the smallest so 
far - inside the nanochannels of porous zeolite AlPOA—h 
single crystals and the authors are considered 0.4 nm di- 
ameter NTs to be ideal one-dimensional quantum hol- 
low wires. The electronic and magnetic properties of 
these ultra-small NTs seem to be drastically different 
from those of large sized NTs. Really given tubes were 
shown to exhibit unusual novel phenomena like diamag- 
netism and superconductivity at low temperatures [23[ 
in addition to several other optical properties and dop- 
ing induced effects. It seems to be the consequence of 
one-dimensionality of given NTs. The direct confirma- 
tion of given conclusion is recent ESR-studies in [24j . 
where ESR measurements on ultra-small single walled 
carbon nanotubes embedded in a SAPO 5 zeolite ma- 
trix with a main point of attention to potentially oc- 
curring CESR (ESR on the electrons in c-band) signals. 
Instead, only one paramagnetic signal was observed of 
symmetric shape at g = 2.0025 on the CNTs in zeolite 
cages, the ESR signal exhibits a predominantly Gaussian 
character at 8.9GHz(4.2K), with a AB pp = 8.3G and 
g = 2.00251. For the free standing SWCNTs, the ESR 
signal takes a more Lorentzian shape at 8.9GHz(4:.2K) 
with ABpp = 5G and with g = 2.00254. So, from 
the above observed distinct changes in the ESR spec- 
tral properties of NTs with and without cages, it was 
inferred in J24], that the ESR signal indeed stems from 
the carbon NTs presenting inside the nanochannels of 
the zeolite. The observed g-value (2.00251) is inconsis- 
tent with the earlier measured values (g = 2.05 — 2.07) 
[25| of CESR on non embedded single- and multiwalled 
NTs of larger diameter. There was observed in [24| a 
dependence of peak-to-peak ESR signal width in con- 
fined SWCNTOSAPO 5 at 70 if on microwave frequecy 
/ suggesting a inhomogeneous (Gaussian) contribution 
to the line broadening. Least-square linear fit of AB PP 
is ABPP{G) = 7.7 +0.018/ '(GHz). Inhomogeneous con- 
tribution to the line broadening is clear evidence, that 
the so called ultra-small carbon nanotubes, studied in 
[24| are not perfect. Nevetheless they can be really con- 
sidered to be quasi-one-dimensional objects. It means, 
that like to t-PA, Peierls transition has to take place. 
Direct indication to given conclusion is g-value, which 
corresponds the only to deep centers in bandgap in any 
carbon materials [ll| , [38| . On the one hand it is indicat- 
ing of the appearance of bandgap in starting metal tubes 
to be consequence of Peierls metal-semiconductor transi- 
tion. On the other hand g — 2.00254 is coinciding in lim- 
its of accuracy of measurements) with principal g-value 
<7i = 2.00255 of axial g-tensor of SSH tt- -solitons in CZS- 
NTs, incorporated in diamond matrix by (111) diamond 
HEIBM and it is near to the g- value of paramagnetic SSH 
7r-solitons in trans-polyacetylene, equaled to 2.00263 [l3j . 
that indicates on the similarity of the structure of tubes 
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incorporated in diamond and embedded in zeolite. At 
the same time angular dependences of g- value and AB PP 
were not observed in [24j. It seems to be additional ar- 
gument to conclusion on insufficient homogeneity of NTs 
in cited work. Really they have different chiralities (au- 
thors report on the presence of three chiralities (5,0) (4,2) 
and (3,3)). On the other hand the observation of spin- 
glass type of temperature dependence of AB PP seems to 
be the indication, that the SSH 7r-solitons are pinned 
and have random distribution along tube axis. For com- 
parison, the temperature dependence of AH PP in NTs, 
produced by HEIBM is quite different and it indicates, 
that AH PP is not determined by unresolved g-anisotropy 
or unresolved hyperfine interaction [l(| ■ It is interesting, 
that pinned solitons can be considered to be similar to 
chemical radicals, possessing by dangling bonds, which to 
some extent is in agreement with interpretation of PC- 
nature to be dangling bond defects in [24| , however they 
are not usual dangling bond point defects, but they are 
topological defects. Random distribution along tube axis 
means in turn, that individual tube with defined chirality 
is also not sufficiently perfect to observe optical quantum 
coherent effects. It means, that Slater principle with one 
preferential direction will be not applicable and optical 
spectra, in particular Raman spectra, will be determined 
the only by local properties of optical centers, that is, 
they can be similar to some extent to those ones in the 
tubes of larger diameters. It is substantial, that main 
Raman modes in the tubes of larger diameters seem to 
be similar in both perfect and in not sufficiently perfect 
NTs, since even in perfect NTs in given case according 
to above consideration there are two preferential direc- 
tion to be existing and one-dimensional quantum coher- 
ent CZSNT-model [lj will be not applicable for descrip- 
tion of optical properties in given tubes. 



B. Quantum Field Effects in Raman Scattering in 
Carbon Zigzag-Shaped Nanotubes 

Scmiclassical evaluation, obtained above has led to 
conclusion, that 7r-subsystem of quasi-lD CZSNTs will 
be inactive in optical experiments. The comparison with 
optical properties of related carbon chain material - car- 
byne - has allowed to predict, that among possible opti- 
cally active topological defects in cr-subsystem of quasi- 
lD CZSNTs the cr-polaron is expected to be prevailing. 
Therefore intensive line 656.8±0.2 cmT 1 in RS-spectrum, 
presented in Figure 1, can be assigned with AFR mode of 
cr-polaron lattice, produced in CZSNTs studied. Further, 
theoretical model, proposed in [l[ allows to insist, that 
the lines, 1215±1 cm" 1 , 1779. 5±1 cm" 1 and 2022.3±0.5 
CTO -i re p resen t themselves the revival part |2| of Rabi 
wave packet in its frequency representation [l|, which is 
result of strong interaction of external EM-field with cr- 
polaron lattice. 

To make given conclusion to be understandable for the 
readers let us reproduce the main moments, concerning 



theoretial model above cited. Each cr-polaron interacting 
with external EM-field in accordance with experiment in 
,27] can be approximated like to guantum dots in Q by 
two- level qubit. Then the Hamiltonian, proposed in the 
work 2] was generalized. The insufficient for the model 
local field term was omitted. (Local field term seems 
to be playing minor role by description of cr-polarons in 
comparison with quantum dots, since size of quantum 
dots is greatly exceeding the size of cr-polarons). The 
apparatus of hypercomplex n-numbers was used. Hyper- 
complex n-numbers are defined to be elements of com- 
mutative ring, given by (|18[) that is, it is direct sum of 
n fields of complex numbers C, n £ N. It means that 
any hypercomplex n-number z is n-dimensional quantity 
with the components k a £ C, a = 0, n — 1, that is in row 
matrix form z is 

z = [k a , h, fc 2 , k n -i], (19) 
it can be represented also in the form 

ra-l 

z = ^ k a n a , (20) 

where 7r a are basis elements of Z n (and simultaneously 
basis elements of the linear space of n-dimensional lines 
and n-dimensional row matrix). They are 

7T = [l,0,...,0,0],7ri = [0,1,..., 0,0], 

...,7r n _i = [0 s 0,...,0,l]. [ ' 

Basis elements 7r a possess by projection properties 
n-i 

Ka^a — ^aS a /3, ^ TTq = 1, Z^a = k a TT a (22) 

In other words, the set of k a £ C, a = 0, n — 1 is the set 
of eigenvalues of hypercomplex n-number z £ Z n , the set 
of {7r a }, a = 0, n — 1 is eigenbasis of Z„-algebra. 

Then the QED-Hamiltonian, considered to be hyper- 
complex operator n-number, for cr-polaron system of in- 
teracting with EM-field ZSCNTs, consisting of n back- 
bones of i-PA chains, which are connected between them- 
selves in that way, in order to produce rolled up graphene 
sheet in matrix representation is [l[ 

[H] = [Ho] + [U F ] + [H*f] + [Uoo] ■ (23) 

The rotating wave approximation and the single-mode 
approximation of EM-field are used. All the components 
in (|23[) are considered to be hypercomplex operator n- 
numbers and they are the following. [Ha] represents the 
operator of the energy of cr-polaron subsystem in the ab- 
sence of interaction between cr-polarons themselves and 
with EM-field. It is 

n-l 

j — 771 
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where a z mj = \a mj ) {a mj \ - \b mj ) (b mj \ is z-transition op- 
erator between the ground and excited states of m-th 
(T-polaron in j-th chain. In other words, transition oper- 
ator technique is used, which was developed in [39| . The 
second term 



n-l 

[Up] = hojY,a + a[ ei Y 

3=0 



(25) 



is the Hamiltonian of the free EM-field, represented in the 
form of hypercomplex operator n-number. Describtion of 
the interaction of cr-polaron sybsystem with EM-field is 
given by the component of the Hamiltonian ([23)) 



a^a + e- ik ™)[e x ] 



n-l 
j—0 m 

(26) 

where g is the interaction constant. The term [H a a] char- 
acterizes intrachain and inter- chain polaron-polaron in- 
teraction. It is given by the relation 

n— 1 n — 1 

\n aa \ = -^E E 4-;i E M N j 

1=0 j=0 m 

71— 1 71—1 

- fi E E4'iiN' E <°»-ui m j 

£=0 .7=0 m 
n — 1 71— 1 

- ^E E4-,iN' E iw (Wwi N j 

n— 1 n-l 

- h J2 E4 2 iiN' E iw ^p-i-.i [ex]'*. 



where ^[i-j] are the energies, characterizing intrachain 
(I = j) and interchain (I j) polaron-polaron interac- 
tion for the excited (£^) and ground states of j-th 



(27) 



are ground and excited states corre- 



chain, \b mj ) , \a. mj ) 
spondingly of m-th (T-polaron of j-th chain. In (HM)) to 
(|27[) [ei] 3 is j-th power of the circulant matrix [ei], which 
is 



ei 
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(28) 



Hamiltonian in the form like to (|2"T1) at n = 1 is usually 
used for description of tunneling between the states with 
equal energies, in particular, for tunneling between quan- 
tum dot states @. Hamiltonian (|27l) at any n describes 
actually the connection between pairs of the states, which 
satisfy the following condition - the first state in any pair 
results from the second state (and vice versa) by time 
reversal. It is known, that for given states Cooper effect 
takes place. Therefore, the application of Hamiltonian 



like to (|27|) is possible for any pair of the states with 
equal energy, which are symmetric relatively time rever- 
sal. By the way, if one omits last term in Hamiltonian 
(|23[) . it goes into rt-chain generalization of well-known 
Tavis-Cummings Hamiltonian [3j. Then the state vec- 
tor of the " NT+EM-field" system was represented [l[ in 
terms of the eigenstates of isolated polaron and photon 
number states in the folowing matrix form 

D*(*)>] = 

n-l 

EiEE «*(*) i<w> +<i(*) iw>)}m j - 

(29) 



j=0 l m 



Here, \b mj ,l) = \b mj ) ® \a mj ,l) = \a mj ) ® \l), where 
|Z) is the EM-field Fock state with I photons, A 3 At), 
B 3 m are the unknown probability amplitudes. Let us 
pay the attention, that matrix function [|W(£))] is single 
state and the representation in the form of sum over j 
is like to representation of complex-valued wave function 
in the form = *i(t) +i^ 2 (t) of two real-valued 
functions. Nonstationary Schrodingcr equation for hy- 
percomplex matrix function in the interaction 
representation is 

d 



[K(i)][|*(t)>], (30) 
where matrix Hamiltonian of interaction [V(£)] is 



[V(t)\ = exp{-[H F ]t){[H a ] + [n aF ] + 

[n aa ])exp(-l[H F ]t) 
n 



(31) 



It leads to the set of difference-differencial matrix equa- 
tions, which in compact form are 

^ t [y m!l (t)] = {-~GJ [a z }-~[E 2 }- 

igVl + l[a z ] exp(i[a z }(ujt - kma))} [* m ,j(t)]+ ^ 
iK]([*m-i,i(*)] + [* m +i,i(*)])» 

where [a z ] is Pauli z-matrix, [E2] is [2 x 2] unit matrix, 
[£] is block matrix 



K] = 2 (&] + ® t^] + (Kil ~ &]) ® M ( 33 ) 

where [£1], [£2] are [n x n] matrices of coefficients, deter- 
mined by (|27| . that is they are 



7i—l n— 1 

[d] = E^-iiNM^] = E^iiN'' ( 34 ) 

2=0 2=0 

where it is taken into account, that in view of axial sym- 
metry [£j 2 ] d° n °t depend on j. Consequently, we have 
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[£x] = [£2] = fe]- Let us remark, that the relax- 
ation processes were taken into consideration in (|32p . It 
was done by means of substituting instead real value ujq 
the complex values with imaginary part iX. It is also 
suggested, that relaxation time r is independent on the 
chain number and it is determined by the value, which is 
reciprocal to A, that is r = 

The solution of hypercomplex equation (|32[) for the 
state vector was obtained in continuum limit and 

it is given by the expression 

[*'(*,*)] = 



$ (h,0)]exp{it([8 l (h)] - gVT+l[ X })}e ihx dh 



(35) 



where x is hypercomplex axis x = [x, x, x], [$'(x, t)] is 

[^(M^ ^-^^U ^fet)], (36) 

In its turn i)] is continuous limit of functional block 

matrix of discrete variable m, which is given by 



[*m,t(*)] = 



(37) 



consisting of two [n x n] matrices of probability ampli- 
tudes 

71-1 

(38) 

71—1 V 7 

3=0 

which arc determined by relationship (|29[) . Further, ma- 
trix [9(h)] in d35J is 

m] = \{(Mh)] + Mh)])®m} 

where [81(h)] and [82(h)] are 



(39) 



[0 1 (/,)] = [6]{2-a 2 (/,+ -) 2 }, 
Mh)] = [&]{2-c?(h-~f} 



(40) 



Here [£1], [£2] are [n x n] matrices of coefficients, defined 
by (HI]). Matrix [x] in ([35]) is 



[x] = i E n] ® exp(-i[cr z ](w - w )t). 
Matrix elements of [$ (ai, i)] are 



(41) 



n-l 

exp {i exp 



0' (/i, 0) exp expihxx 



lixqji 



(42) 



(^■(ft) - gVl - lKj(h))}dh, 



where @ l q (h, 0),-&j(h), Hj(h) are determined by eigenval- 
ues k Q e C, a = 0, n - 1 of (/i, 0) , 0(/i) and x(/i) , which 
are considered to be n-numbers. They are 



I 1 

0<(M) = -k,($'(/i,0)) - - V<I><(M)exp 
* n n * — ' J 

3=0 



2nqji 



(43) 



1 1 n 9 ' ' 

0j(h) = -kMh)) = -^^)^P— , (44) 



r=0 



n-l 



= ^MxW) = I^^-^exp^. (45) 



!-=0 



Then the hypercomplex solution can be represented in 
the form of sum of n solutions for n chains, that is, hy- 
percomplex n- number &(x,t) is 



n-i 

$ l (x,t) = '£¥ q (x ) t), 

q=0 

where the solution for q-th chain (x, t) is 



§ I ,(* J t) = £^ p (a: ) i)[e 1 ]'\ 

p=0 



(46) 



(47) 



in which the matrix elements & qp (x,t) are determined 
by (|42l) . The relationship (|47|) by taking into account 
(1551 - (HHJ) determines Rabi-wave packet, which prop- 
agates along individual chain of zigzag NT. It is now 
clear, that subsequent analysis of Rabi-wave packet dy- 
namics for individual NT-component will be the same 
(by rescaling of parameters), that in [2]. The parame- 
ters can in principle be obtained by detailed comparison 
of results obtained with experiments on optical absorp- 
tion, reflection or Raman scattering. It is the subject 
for futher work. For qualitative analysis we can use 
arbitrary parameters of the task. Then from equation 
(1471) follows temporal dependence of the integral inver- 
sion for j-th chain in CZSNT and spectral dependence 
of the RS-signal amplitude, which explain the appear- 
ance of a number lines, which are additional to AFSWR 
lines in the spectra observed. Really temporal depen- 
dence of the integral inversion gives in implicit form the 
way for comparison of theoretical results with any sta- 
tionary optical experiments in QD-chain or like them, 
including quasi-lD CZSNTs, with aforesaid initial state. 
It is sufficient to make a Fourier transform of given tem- 
poral dependence. It will be proportional to signal am- 
plitudes of infrared (IR) absorbtion, IR-transmittance, 
IR-rcflection or Raman scattering, since they are deter- 
mined by population difference. It means, that dynam- 
ical nonstationary properties of optical systems can be- 
come apparent by conventional stationary registration of 
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the spectra. It is very similar to the well known in sta- 
tionary ESR-spectroscopy situation, where the signals of 
the centers with long relaxation times can be registered in 
nonstationary regime, which can lead to zeroth absorp- 
tion (or appearance of the signals in inverse phase like 
to maser effects), that was practically used in very many 
ESR-studies to unravel complicated overlapped spectra 
(however without theoretical explanation of given effect). 
At the same time there is difference of principle between 
classically considered nonstationary regime with Rabi os- 
cillations by classical EM-field consideration and between 
QED-consideration. If EM-field is considered to be clas- 
sical field, then operating point by the registration of the 
centers with long relaxation times can move along Rabi 
oscillation curve, leading to radical changes in amplitude 
of the signals, including the appearance of irradiation in- 
stead absorption. Naturally in given process the energy 
of EM-field source is partly comes back from spin sys- 
tems in the same quality (in EM-field form) instead its 
transformation in phonon energy. Given process can be 
accompanying by enhanced noice level, however without 
appearance of any additional lines. QED-consideration 
shows, that the process of Rabi waves' formation is de- 
termined essentially by interaction of optical centers with 
photon subsystem, at that stationary regime is achieved 
fastly, giving usual stationary optical signals, however af- 
ter a time the oscillations emerge again, and given revival 
part leads to additional lines in high energy part of spec- 
tral dependence. Therefore, the absorption (scattering) 
process by presence of Rabi waves is always nonequilib- 
rium process in the whole system [EM-field + matter] 
and it can be realized even at relatively short relaxation 
times of optical centers, if constant g of the interaction 
of optical centers with EM-field is large. 

Thus, the CZSNT-QED-model predicts the appear- 
ance of additional lines, for which Rabi waves are respon- 
sible, in stationary optical spectra of quasi-lD CZSNTs. 
It is understandable, that Rabi waves can be registered 
in any ID-systems by stationary optical measurements, 
if the electron-photon coupling is rather strong. 

The identification of the lines 1215±1 cm -1 , 1779. 5±1 
cm -1 , 2022. 3±0. 5 cm -1 (and corresponding lines in 
boron ion modified sample) with Fourier-image of revival 
part of Rabipacket is confirmed by the following. It is 
well known [8j , that in the case of point absorbing cen- 
ters Rabi frequency is linear function of the amplitude 
of oscillating EM-field. I was experimentally confirmed 
in [4lj] , (42[. Given dependency takes also place for the 
center of Rabi wave packet, that follows from the anal- 
ysis of Fourier transform of temporal dependence of the 
integral inversion. Further, it is evident, that amplitude 
of laser wave, penetrating in IBMR, is lesser by excita- 
tion from unimplanted side. It is consequence of some 
absorption in the unimplanted volume of diamond crys- 
tal. We see, that really, the values 1569 ±3 cm -1 , 1757±5 
cm" 1 of two high frequency components are substantially 
less in given case, than 1779. 5±1 cm -1 and 2022. 3±0. 5 
cm -1 , observed by excitation from implanted side, Fig- 
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FIG. 4: Dependence of 2D-band position in graphene on ex- 
citation energy, presented in [i^] (Figure 5) 



ure 2, Figure 1. Relative frequency changes are 1,151 and 
1.134 ±0.003 for the pairs [2022.3±0.5 cm -1 , 1757±5] 
and [1779. 5±1 cm" 1 , 1569 ±3] correspondingly and it is 
seen, that they are close to each other, at that larger fre- 
quency undergoes larger change in correspondence with 
theoretical analysis. 

The lines at 354.6, 977.1 (±1 cm -1 seem to be as- 
signed with two AFESWR modes. The substantial de- 
crease of relative intensity of 641.8 ±1 cm -1 mode in 
comparison with 656.8±0.2 cm -1 mode testifies in favour 
of given assignment. It is seen, that AFESWR-splitting 
is rather large and it has thr same order of values with the 
splitting between two polaron vibronic levels. It means, 
that linear AFESWR-theory (27|, which predicts a set of 
equidistant AFESWR-modes, arranged the left and the 
right of central mode, can be used the only to obtain 
approximately the average value of AFESWR-splitting. 
Really, in given case AFESWR-modes are not equidis- 
tant, they are shifted on different distances 335.3, 287,2 
cm" 1 from main AFESWR- mode. At the same time av- 
erage value of AFESWR-splitting is 311.3cm -1 and it is 
close to the value in 300 cm -1 expected in accordance 
with linear AFESWR-theory by taking into account ex- 
perimental AFESWR-splitting value in carbynes. 

Therefore, we obtain the direct proof of assigment of 
the lines 656.8±0.2 cm -1 and 641.8 ±1 cm" 1 with a- 
polaron lattice. 

To explain differences by the change of the direction of 
the excitation wave propagation, we have to take into 
consideration the following. The possibility of SWR- 
excitation in ID-systems is strongly dependent on the 
geometry of experiment, detrmined by directions of chain 
axis, vectors of intracrystalline and external magnetic 
or electric fields (by ferromagnetic or ferroelectric SWR 
study correspondingly). It is confirmed by ferromagnetic 
SWR study in carbynes, [26[. The axes in CZSNTs are 
not linear in the end of ion run and they are genera- 
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FIG. 5: Dependence of 2D-band position in graphene on elec- 
trical component of oscillating excitation field 



trixes of the figure of onion-like shape, that provides for 
necessary geometry for AFESWR-excitation of er-polaron 
lattice. Moreover the appearance of very broad line, Fig- 
ure 2, seems to be indication on the excitation of the 
Frohlich movement of cr-polaron lattice itself. The pres- 
ence of Frohlich sliding of cr-polaron lattice allows to ex- 
plain qualitatively the appearance of "hysteresis" in spec- 
tral dependences relatively changing of the direction of 
excitation wave propagation into opposite direction. It 
follows from energy law conservation position. Really, 
moving cr-polaron lattice possesses by kinetic energy. It 
means, that it is required lesser energy value to excite 
the local polaron vibration mode in correspondence with 
observation. 

Therefore CZSNTs incorporated in diamond matrix 
represent themselves the example of the system, which 
strongly interact with EM-ficld. Experimental detection 
of Rabi wave packets confirms on the one hand the the- 
ory, developed in [l| and it is also the first experimental 
confirmation of Rabi wave phenomenon, predicted in Q , 
on the other hand. It means also, that semiclassical de- 
scription of spectroscopic transitions in CZSNTs and in 
the systems like them cannot be appropriate. It substan- 
tially raises the practical concernment of QED-theory. 

The requirement of perfect one-dimemensionality is 
substantial for applicability of Slater principle to CZS- 
NTs, at the same time the Rabi wave phenomenon seems 
to be general and can be observed in 2D and 3D-systems. 
For instance, the first-order RS-peak in the vicinity of 
1580 cmT 1 and so called second-order lines around 2480, 
2700 (2D band), 3250 cm" 1 by X exc = 488 nm were ob- 
served in graphene 40] , at that the intensity of second- 
order 2D-band was found to be substantially exceeding 
the intensty of the corresponding first-order line. More- 
over theoretical dependence of 2D-band position on ex- 
citation energy, presented in [40j by Figure 5, is far from 



experimental dependence, see Figure 4, which is repro- 
duction of Figure 5 from |4(| . At the same time the sug- 
gestion, that the second-order lines are spectral mapping 
of revival part of Rabi wave packet, corresponding to the 
first-order line is agreeing with experimental data in |40| 
very well, see Figure 5 in our paper. It wasobtained from 
experimental data in [40| by choosing of ^-coordinate in- 
stead the energy the electrical component of oscillating 
excitation field. It is seen, that given experimental de- 
pendence is linear in full correspondence with expected 
for the center of Rabi wave packets. It is analogue of well 
known in radiospectroscopy linear dependece of Rabi os- 
cillation frequency on magnetic component of oscillating 
excitation field, see for instance [4l[, [42| . 

Therefore given result is experimental evidence of the 
formation and propagation of Rabi wave packets in 
graphene, being to be the example of 2D-system. Simi- 
lar explanation can be proposed for transitions in heavily 
boron-doped diamond in the region of 1400-2800 cm -1 
[43l | . It means, that really the theory of Rabi wave, pre- 
sented in for ID-systems can be developed and gen- 
eralized for 2D- and 3D-systems. Consequently, it seems 
to be evident, that Rabi wave packets can also be iden- 
tified in free standing NTs by the production technology 
improvement. 



V. CONCLUSIONS 

Thus, QED-model for multichain coupled qubit quasi- 
1D system, proposed in [l[ , is confirmed by Raman scat- 
tering studies of carbon zigzag-shaped nanotubes, pro- 
duced by HEIBM of natural diamond single crystals. 
New quantum optics phenomenon - Rabi wave packet 
formation, predicted in Q , has been experimentally iden- 
tified for the first time. It is shown, that Raman spectra 
in quasi-lD CZSNTs are quite different in comparison 
with well known Raman spectra in 2D those ones. They 
characterized by semiclassical consideration by the only 
one vibronic mode of Su-Schriffer-Heeger cr-polaron lat- 
tice with peak positions 656.8±0.2 cm -1 , for instance, in 
copper implanted sample (by excitation from impantcd 
side of the sample) instead of longitudinal and transverse 
optical phonon G + and G~ modes and the out-of-plane 
radial breathing mode, which are observed in Raman 
spectra of 2D single wall nanotubes. It is consequence 
of 2D - ID transition in all physical properties of nan- 
otubes. 

It is shown, that strong electron-photon coupling takes 
place in quasi-lD CZSNTs by interaction with EM-field 
and quantum nature of EM-field has to be taken into 
account. It has been done for the first time in station- 
ary spectroscopy at all. All optical spectra, in particular, 
Raman spectra are registered by usual stationary mea- 
surement technique in nonequilibrium conditions, which 
are the consequence of Rabi wave packets' formation. 
It leads in its turn to appearance of additional lines 
with peak positions 1215±1 cm -1 , 1779. 5±1 cmT 1 and 
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2022. 3±0. 5 cm" 1 in copper implanted sample and addi- 
tional lines with peak positions at 1212. 3±1, 1772. 5±1, 
2011±0.5cto~ 1 in boron implanted sample (by excita- 
tion from implanted side in both the samples), which are 
corresponding to revival part of inversion dependence of 
joint EM-field + matter system in the frequency repre- 
sentation. 

The substantial qualitative and quantitative depen- 
dence of Raman spectra on changing of the direction of 
excitation wave propagation into opposite direction has 
been revealed. Instead of the RS-lines with peak posi- 
tions at 656.8±0.2 cm" 1 , 1215±1 cm" 1 , 1779. 5±1 cm" 1 
and 2022. 3±0. 5 cm -1 , which are obseved by excitation 
from implated side of the copper imlanted sample, RS- 
lines with peak positions at 1328.7 cm -1 , 354.6, 641.8, 
977.1 (±1 cm -1 ), 1274.1 ±2 cm" 1 and more weak pro- 
nounced lines at 1569 ±3 cm" 1 , 1757±5 cm -1 were ob- 
served. Moreover, all the lines listed are superimposed 
with very broad (its linewidth is 1720±20 cm" 1 ) asym- 
metric line with peak position at 1160±10cra . RS-lines 
with peak positions at 354.6, 641.8, 977.1 (±1 cm" 1 ) 
were identified with AFESWR-modes. RS-line with peak 
position at 1274.1 ±2 cm" 1 and the lines at 1569 ±3 



cm" 1 , 1757±5 cm" 1 were reffcred to revival part of in- 
version dependence of joint EM-field + matter system 
in the frequency representation. More low frequency val- 
ues in comparison with those ones observed by excitation 
from implanted side of the sample are connected with de- 
crease of exciting light intensity to be the consequence of 
additional absorption in unimplanted region of the sam- 
ple. RS-line with peak position at 1328.7 cm" 1 corre- 
sponds to optical phonon in diamond lattice, its value 
indicates on some softening of optical phonon mode, tak- 
ing place by copper implantation. It is interesting, that 
softening is absent by boron implantation. It seems to 
be display of size effect relatively effective radius of im- 
plated ions. Very broad asymmetric line with peak po- 
sition at 1160±10 cm" 1 is attributed to Frohlich sliding 
of cr-polaron lattice. The appearance of sliding allows to 
explain qualitatively on the basis of energy conservation 
law the effect of "hysteresis" in the frequency value of 
main AFR mode by changing of the direction of excita- 
tion wave propagation into opposite direction resulting 
in its decrease to 641.8 ±1 cm" 1 ) instead of 656.8±0.2 
cm" 1 . 
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